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ABSTRACT 

The soft gamma repeater SGR 1806—20 is most famous for its giant flare from 2004, 
which yielded the highest 7-ray flux ever observed on Earth. The flare emphasized 
the importance of determining the distance to the SGR, thus revealing the flare's 
energy output, with implications on SGRs energy budget and giant flare rates. We 
analyze x-ray scattering echoes observed by Swift /XRT following the 2006 August 6 
intermediate burst of SGR 1806—20. Assuming positions and opacities of the molecular 
clouds along the line-of-sight from previous works, we derive direct constrains on the 
distance to SGR 1806—20, setting a lower limit of 9.4 kpc and an upper limit of 18.6 
kpc (90% confidence), compared with a 6—15 kpc distance range by previous works. 
This distance range matches an energy output of sa 10^^ ergs~^ for the 2004 giant 
flare. We further use, for the first time, the x-ray echoes in order to study the dust 
properties in molecular clouds. Analyzing the temporal evolution of the observed flux 
using a dust scattering model, which assumes a power-law size distribution of the dust 
grains, we find a power-law index of — 3.3]to 7 (I''') ^-nd a lower limit of O.l/xm {2a) 
on the dust maximal grain size, both conforming to measured dust properties in the 
diffused interstellar medium (ISM). We advocate future burst follow-up observations 
with Swift, Chandra and the planned NuSTAR telescopes, as means of obtaining much 
superior results from such an analysis. 
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1 INTRODUCTION 

Soft Gamma Repeaters (SGRs) are objects emitting soft 
gamma-ray and hard x-ray bursts at irregular intervals, 
as well as a persistent x-ray emission (for a review see 
Woods & Thompson 2006 and references therein). Bursts 
are typically short (~ 0.1s) and are gathered within active 
periods that last between a few weeks to several months, 
followed by years of quiescence (Kouveliotou 1998). Bursts 
are commonly classified according to their peak luminos- 
ity, from the most common flares reaching io*^~*^ ergs^^, 
up to the rare 'giant' flares reaching ~ 10*'' ergs~^. SGRs 
are believed to be magnetars — neutron stars with a surface 
magnetic held of ~ 10^^ G, which serves as the energy source 
of the bursts and the persistent emission (Paczynski 1992; 
Duncan & Thompson 1992; Kouveliotou et al. 1998). 

SGR 1806—20, lying in the direction of the Galactic 
center behind a veil of 15 — 30 magnitudes of optical ex- 
tinction (Corbel & Eikenberry 2004, hereafter CE04), is one 
of a handful known SGRs. Its resume includes the flrst 
SGR to be observed, on January 7 1979 (originally classi- 



fled as a gamma ray burst, Mazets & Golenetskii 1981) and 
the subject of the flrst SGR spindown rate measurement 
(Kouveliotou et al. 1998) — a major milestone in the accep- 
tance of the magnetar hypothesis. Yet it is most famous for 
producing the most energetic giant flare observed to date; 
on 2004 December 27 it emitted a 0.1s flare with an esti- 
mated total energy of 2 — 5 x 10*'' (^^s^pc) ^^^'^^ 
the distance to the SGR (Palmer et al. 2005; Hurley et al. 
2005). The emitted energy in this event is evaluated 100 
times higher than the next most energetic SGR recorded 
event (assuming d = 15 kpc). 

An energy of i5 ~ 10*® erg is at odds with some of 
the observations. First, a naive rate of yr~^ flares with 
similar energy per SGR^ is ruled out by failure to observe 
corresponding population of extragalactic SGRs (detectable 
to ~ 30 Mpc, Palmer et al. 2005; Nakar et al. 2006). More- 



^ based on a single E 10*'' ergs event within 30 years of ob- 
serving ~ 5 SGRs, as of 2004 
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over, even a rate of a single giant flare per SGR life span, 
within 2a of the observed rate obtained by a careful statis- 
tical treatment, has only marginal agreement with observed 
extragalactic rates (Ofek 2007). Second, the surface mag- 
netic fleld corresponding to observed properties of SGRs 
(e.g. spindown rate) is B « lO^^G, matching an external 
magnetic energy of ~ lO'^^erg, comparable to the energy out- 
put of a single giant flare. Since the source of both the energy 
reservoir powering giant flares and the reservoir responsi- 
ble for the persistent emission is thought to be the surface 
magnetic field of the SGR, one might expect a significant 
transformation in the observed spectral and temporal emis- 
sion parameters of SGR 1806—20 following the giant fiare. 
This expectation is not met by observations'^ (Woods et al. 
2007; Esposito et al. 2007). A possible explanation for both 
the rates discrepancies and the unchanged emission features 
is a distance shorter than the commonly assumed 15 kpc to 
SGR 1806—20, matching a less energetic flare output. 

Employing different approaches, several papers from re- 
cent years suggest distance ranges within 6—15 kpc to SGR 
1806—20. The emerging factor of nearly 3 in the SGR's dis- 
tance estimates translates to nearly an order of magnitude 
difference in its emitted energy. CE04 used CO emission lines 
and NH3 absorption features from molecular clouds along 
the line of sight to find the clouds radial velocities, inferring 
two possible locations per cloud, one in front of the Galactic 
center and one behind it. Accounting for the optical extinc- 
tion of the star powering the nebula LBV 1806—20 they 
determined a distance of 15.llit'3 kpc {2a) to the cluster 
containing LBV 1806-20. They associated SGR 1806-20 
with this cluster due to its angular proximity of 12" to LBV 
1806-20 and the match between SGR 1806-20 x-ray ab- 
sorption and the IR extinction towards the cluster mem- 
bers. Figer et al. (2004) measured radial velocities using ab- 
sorption lines from LBV 1806—20 and nearby stars, which 
translated to a distance of 11.8 kpc. Bibby et al. (2008) 
spectroscopically classified several stars which were identi- 
fied as members of the cluster of LBV 1806-20 by CE04 
and Figer et al. (2005). Based on their absolute magnitude 
calibration and near-IR photometry, as well as isochrones 
fit to the cluster's age, they obtained a cluster distance of 
8.7l}'5 kpc. As opposed to the above associative distance 
estimations, Cameron et al. (2005) gave a more direct es- 
timate. They identified the decaying bright radio afterglow 
of the December 2004 giant flare a week after the burst. 
Based on absorption features of intervening interstellar neu- 
tral hydrogen clouds, they constrained SGR 1806—20 dis- 
tance to within 6.4—9.8 kpc. McClure-Griffiths & Gaensler 
(2005) accepted the lower limit of ~ 6 kpc but rejected their 
upper limit, disqualifying the association of the absorption 
feature used to set this limit with SGR 1806—20. 

We present a new direct estimate of the distance to SGR 
1806-20 based on dust scattered x-ray observations, from 
which we also extract properties of the dust along the line 
of sight. The scattering of x-rays from dust grains in the ISM 

^ These works do report changes in temporal characteristics of 
the SGR following the giant flare, suggesting reorganization of 
the magnetic field. Similar reports followed the giant flare of SGR 
1900+14, e.g. G6gii§ et al. (2004). However, the loss of a signifi- 
cant portion of the total energy reservoir seems to justify a more 
dramatic change. 



was first considered by Overbeck (1965). For an x-ray source 
with a varying intensity (e.g. a short burst), if the dust spa- 
tial distribution is known, one can use the time delay be- 
tween the direct signal and the scattered signal to constrain 
the distance to the x-ray source (Trumper & Schonfelder 
1973). This was first applied for constraining the distance 
to the X-ray binary Cyg X-3 (Predehl et al. 2000). 

Analysis of x-ray halos around Galactic and extra- 
galactic sources have been used to constrain the proper- 
ties of dust grains in the ISM (Mauche & Gorenstein 1986; 
Predehl & Schmitt 1995), with results conforming to the 
dust model by Mathis et al. (1977). Such analyses exploit 
the dependence of the scattering cross section on the grain 
properties as well as the dependence of the halo radial pro- 
file on the positions of the scattering grains along the line of 
sight. 

In the case of both a short burst and a thin dust scat- 
terer, the halo is replaced by a ring which radius increases 
with time. Measuring the expansion rate of such observed 
rings from GRBs has been used to derive distances to Galac- 
tic dust clouds (Vianello et al. 2007) as well as the time of 
the original burst (Feng & Fox 2010). Tiengo et al. (2010) 
combined observations by both Swift fKKT and XMM- 
Newton/EPIC of rings following bursts of the anomalous 
X-ray pulsar 1E1547. 0-5408, and used several different mod- 
els for the dust composition and grain size distribution to 
fit the intensity decay of each ring as a function of time 
and energy, in order to obtain constrains on the distance 
to the X-ray source. They concluded that in the absence of 
independent constraints on the distance of the source or the 
scatterer, their analysis is highly sensitive to the size of the 
largest dust grains, as reflected by each of the models used. A 
dust-scattered halo surrounding an SGR was flrst reported 
by Kouveliotou et al. (2001) for SGR 1900+14, with data 
quality insufficient for further analysis. 

The dust along the line of sight to SGR 1806-20 is 
concentrated within molecular clouds (e.g. CE04). As op- 
posed to the diffuse ISM, where dust properties are ef- 
ficiently probed using UV and longer wavelengths, the 
dust in cores of molecular clouds is not easily accessible 
due to the clouds high optical depth at these wavelength. 
Moreover, observational evidence (Carrasco et al. 1973; Jura 
1980; Goldsmith et al. 1997; Stepnik et al. 2003; Chiar et al. 
2007; Winston et al. 2007; Schnee et al. 2008; Butler & Tan 
2009) and theoretical considerations (e.g. Ossenkopf 1993; 
Weidenschilling & Ruzmaikina 1994; Ormel et al. 2009) in- 
dicate that the high density in molecular clouds leads to 
grain coagulation that alters the dust grain size distribution 
from the one observed in the diffuse ISM. Since molecular 
clouds are optically thin to hard x-rays, small angle x-ray 
scattering provides a unique tool to directly probe grain 
size distribution at their cores. Despite this virtue, x-rays 
have not been used yet to probe dust properties in molecu- 
lar clouds. Here we seize the opportunity to employ, for the 
first time, x-ray echoes for this purpose. 

On 2006 August 6, an intermediate burst of SGR 
1806—20 was observed (Hurley et al. 2006b, a). Anticipating 
delayed dust scattered echoes. Swift fKKT took a dozen ob- 
servations during the two weeks following the burst. A pre- 
liminary analysis of the first two observations by Goad et al. 
(2006) reported an expanding halo due to dust scattering. 
We reanalyze Swift observations, finding the halo fiux profile 
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of each observation. Assuming properties of molecular clouds 
along the line of sight to SGR 1806-20 as reported in CE04, 
we use the observation profiles to constrain the distance to 
SGR 1806-20. We then assume, instead of the CE04 dust 
distribution, a single predominant dust screen along the line 
of sight and a power-law distribution for the dust grain size, 
and use the observations to constrain the grain size distribu- 
tion within the intervening molecular cloud that dominates 
the scattering. 

This paper is organized as follows. Observations are de- 
scribed in Section 2. Section 3 reviews the dust scattering 
model we use in our analyses, which are presented in Section 
4. In section 4.1 we assume the distance of the scattering 
clouds is known from CE04 and constrain the distance to 
SGR 1806—20. In section 4.2 we relax our assumptions re- 
garding the distance to the scattering clouds and study the 
properties of the scattering dust. We discuss implications of 
our analyses for future observations in section 5 and draw 
conclusions in section 6. 



2 OBSERVATIONS 

SGR 1806—20 intermediate burst of 2006 August 6 was first 
reported by Hurley et al. (2006b, a). It comprised six sepa- 
rate bursts over «120 seconds, with the dominant one last- 
ing for ~30 seconds and characterized by a measured flu- 
ence of 2.4 x 10^* ergcm^^ and an optically thin thermal 
bremsstrahlung (OTTB) spectrum of kT = 20 keV within 
the Konus-Wind range of 20—200 keV (Golenetskii et al. 
2006). Following this report. Swift /XRT took 12 observa- 
tions in Photon Count mode, starting 30 hours after the 
burst and ending 14 days later, corresponding to observation 
ID 00035315002-00035315013. For our analysis we used the 
four earliest observations (see Table 1, for 1st and 2nd obser- 
vations see Figure 1). The other observations did not show 
a significant signal and were only used implicitly for back- 
ground consistency check. Swift /XKT is sensitive to photons 
in the energy range 0.2—10 keV, and the dust scattered sig- 
nal in the observations is only evident in the range 2.5—8.5 
keV. Below this range the optical depth is larger than unity 
and the assumption of a single scattering does not hold. In 
order to allow energy dependent analysis without diluting 
the signal too much we split each observation data into two 
energy bands: 2.5—4.5 and 4.5—8.5 keV. We created expo- 
sure maps for each observation and each energy band with 
the FTOOLS task XRTEXPOMAP, to correct for the vi- 
gnetting, the dead detector areas, and the excluded regions. 
For each observation we extracted a radial profile centered at 
the position of SGR 1806-20 given by FTOOLS task XRT- 
CENTROID. Radial bin width of 8 pixels (19") was chosen 
so it overlaps with average Swift /XICT point-spread function 
(PSF) of ^ 20" (half power diameter. Burrows et al. 2005). 
We estimated the background photon count per energy band 
for each observation by taking the average photon count at 
a radial distance where no dust scattered signal is expected 
and adjusting for the area of each radial bin and for its av- 
erage exposure time derived from the exposure map. The 
background thus measured agrees across all 12 observations 
with a photon count C = 1.40 ± 0.13 x 10~^ phs~^ for the 
2.5-4.5 keV band and C = 9.57± 1.20x10"^ phs"^ for the 
4.5—8.5 keV band over the XRT detection area. 



The profile cannot be explained by the scattered signal 
and the background alone, and is found to consist of a third 
component — a fading halo around the SGR. This halo is 
wider than the central source PSF and therefore also con- 
tains one or more of the following contributions: an integra- 
tion over dust scattered rings originated from fading post- 
burst source emission at times between the burst and the 
observation; dust scattered rings unresolved from the source 
due to dust clouds adjacent to the source; and multiple scat- 
tering. Lacking knowledge about the post-burst source lu- 
minosity as a function of time and confined by the XRT 
resolution, we model this halo as a two parameters King 

function C* ( 1 + (;^) I ■ This component has no signif- 
icant effect on our background evaluations since the King 
function is negligible at the radial distance used for back- 
ground measurement. 



3 DUST SCATTERED X-RAY RINGS 

An SGR flare can be considered a special case of a varying 
source, where the flare is approximately a delta function in 
time, and the result of an intervening dust screen is a ring 
expanding with time. For small scattering angles and a single 
scattering screen, geometry dictates 



where 0{t) is the measured ring angle with respect to the line 
of sight, d is the source distance, x is the screen's location 
expressed as a fraction of d, and t is the time passed since 
the direct flare observation. 

We fit the observations to a model based on 
Smith & Dwek (1998), with differential scattering cross 
section as described in Rivera-Ingraham & van Kerkwijk 
(2010). The original model calculates the halo observed by 
the scattering of x-rays from a persistent source over a con- 
tinuous distribution of intervening dust. We are interested 
in the rings produced by the scattering of an impulsive emis- 
sion over one or more molecular clouds, namely a discrete 
distribution of intervening dust. Our configuration is simpler 
and it allows us to derive an analytic solution. We therefore 
modify the original derivation by Smith & Dwek (1998) to 
obtain an analytic expression for the fiux in a ring produced 
by the scattering of a flare on a single thin dust screen, i.e. 
the scattered flux observed at an angular distance Sobs(i) 
from the central source, as a function of photon energy. 
The modification is straight forward and we therefore do 
not repeat the derivation of Smith & Dwek (1998) here, and 
just highlight the modifications. Keeping the original nota- 
tions, we replace the dust distribution along the line of sight 
with the dust column density across the scattering screen, 
A'dust and the dust grain size distribution with A^a, such that 
A^dust = / Na.da where a is the radius of the dust grain. We 
further replace the fiux reaching the scatterer, used in the 
original derivation, by the flare's total fluence reaching it. In 
addition, the solid angle which the scattered observed sig- 
nal occupies, dO! , is 2Ti9{t)d9 where dO is the change in the 
observed angle during the observation time dt. Using Equa- 
tion 1 we find dO.' = g{x)dt where g {x) = The 
observed scattered fiux per photon energy E, Fe (energy 
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Table 1. Swift /XRT observations used in the analysis. Times are given in hours, start and end times are relative to the burst trigger. 
Total duration is the time elapsed between the start and end time while exposure time is the net time the detector was actively collecting 
signal. 



Obs. ID 


Start time 
[hr] 


End time 
[hr] 


Total duration 
[hr] 


Exposure time 
[hr[ 


00035315002 


30.50 


31.17 


0.66 


0.66 


00035315003 


82.22 


104.81 


22.59 


3.89 


00035315004 


112.94 


129.17 


16.23 


1.89 


00035315005 


141.93 


153.29 


11.36 


1.88 



08/07/2006, 0.66 hr expo: 



-400 - 
-600 - 



-800 -600 -400 -200 200 400 600 800 
arcsec 



Figure 1. The dust rings seen in the 1st (left, «31 hr after the burst) and 2nd (right, f»90 hr after the burst) Swi/t/XRT observations, 
(energy range 2.5—8.5 keV). 



observed in a ring per photons of energy E per unit area 
per unit time dt) is thus: 



Fe i9obs[t]) = Se 



(1 



7,r I da 

""^[dn 



da 



(2) 



where Se is the flare's direct (vs. scattered) observed fluence 
per photon energy E and da{E,a,9sca.t)/dQ is the differen- 
tial scattering cross section of an x-ray with energy E hy a 
dust grain with radius a at an angle Sscat. 

As scattering angles are small, we substitute 6'scat ~ 
9obs/ {1 — x). For the dust grain size distribution we assume a 
power-law A^a = Aa'' where amin ^ a ^ ctmax (Mathis et al. 
1977). For the differential cross section we use the Gaus- 
sian approximation of the Rayleigh-Gans theory, valid for 
small scattering angles and energies above IkeV, following 
van de Hulst (1957): 



da{E,a,6s, 
dQ 



Cdustts exp 



a2 

^scat 
2^scat 



(3) 



where Cdust depends on dust components atomic charge, 
mass number, density, and scattering factor, while Oscat is 
given by Mauche & Gorenstein (1986) 

1 



{a,E) = 10.4 



(£'/lkeV)(a/0.1^im) 



(4) 



Inserting the cross section and dust grain size distribution 
we obtain 



-pE(^obs, 9, fimax, Cimin) — 

B,/5;;-;a'«exp|-i(^)'(^ 



i/jm 



E 



(5) 



da 



where a ; 



and 



Be = SECdustA—p{l 



-,9+6 



/X 



(6) 



The differential cross section (Equation 3) implies that a 
dust grain of size a scatters effectively below a characteristic 
angle 9 (a), above which the differential cross section decays 
rapidly. Thus, for a grain size distribution more gradual than 
a~^, regardless of its exact form, the scattering at an angle 
Sscat is dominated by the grain size matching Equation 4, 



^Iscat (6scat, E) 



0.1 



(E-ZlkeV) (6»scat/10.4arcmin) 



/xm 



(7) 



The observed scattered flux can then be approximated as 
Fe{9) oc Aa g^j^t a Jcat • If the grain size distribution is lim- 
ited to the range amin ^ a ^ ^max then at small angles 
(early time) a scat > flmax and the scattering is dominated 
by the largest grains, namely it is constant in angle (time) 
and Fe oc Na^^^^a"^^^. At large angles (late time) when 
a scat < dmin the Scattering is dominated by the Gaussian 
tail of the smallest grains cross-section, and the flux falls 
exponentially. Therefore, defining 9 — (1 — a;) 6'scat, Equa- 
tion 5 can be approximated as 



Fe{9) 



const 

t^obs 

exp 



^obs ^ 9 (Oniax) 

9 (amax) < fi'obs < 9 (amin) (8) 
[- ^obs > 9 (amin) 



It is useful to express this approximation in terms of the flux 
dependence on time, using Equation 1: 
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i const t <^t (amax) 

f(a„,ax) «t«f(a„,in) (9) 

where t {a, E) = 2c{i-x) These expressions describe 

three flux regimes — (i) a constant set by scattering over 
the largest grains, (ii) A power-law decay where the size of 
the grains that dominate the scattering vary with time and 
(iii) an exponential decay set by the scattering cross section 
of the smallest grains. Therefore, observing the evolution 
of Fe with time is a direct probe of the dust grain size 
distribution. Omax is probed by observations covering the 
time of transition between the constant and the power-law 
regime, while the index q is obtained by power-law regime 
observations. Probing Omin requires a burst bright enough 
for the scattered signal to overcome the background at the 
low end of the power-law regime. 

We note that Equation 10 in Draine (2003), which is an 
analytical approximation for the differential scattering cross 
section describing his model of x-ray scattering by dust, can 
be closely matched by the above scattering model by choos- 
ing Omax = 0.33 /im and q = — 3, which are typical values 
for ISM dust. 



Table 2. Molecular clouds along the line of sight to LBV 
1806-20, taken from CE04. The bottom three clouds, lack- 
ing unambiguous distance determinations, are quoted with both 
near/far distances (Icr errors). CE04 do not determine whether 
the three last clouds in this table are in front or behind the Galac- 
tic center, allowing two distances to each cloud. In the text we 
refer to all three clouds being in front of the Galactic center as 
the 'near' configuration, while all three clouds being at their far 
distance is the 'far' configuration. 



Name 


Av 


Distance 




[magnitude] 


[kpc] 


MC-16 


8.6±1.7 


4.5±0.3 


MC4 


3.9±0.8 


0.2±0.5 


MC13A 


10.0±2.0 


15.1±0.9 


MC13B 


3.0±0.6 


4.5±0.3 


MC24 


5.6±1.1 


3.0±0.5 


MC30 


0.8±0.2 


3.5±0.35 


MC38 


1.1±0.2 


4.2±0.3 


MC44 


1.6±0.3 


4.5±0.3 


MC73 


6.0±1.2 


5.7/11.0±0.15 


MC87 


1.4±0.3 


6.1/10.6±0.15 


MC94 


0.5±0.1 


6.2/10.5±0.15 



4 ANALYSIS & RESULTS 

To obtain our two objectives — constraining the distance to 
SGR 1806—20 and the properties of the dust in the inter- 
vening molecular clouds, we take two different approaches 
in analyzing the XRT observations. First (Section 4.1), we 
assume positions and visual extinctions of molecular clouds 
along the line of sight based on CE04 in order to get an 
estimate of the distance to SGR 1806—20. We test the de- 
pendence of this analysis on a dust scattering model and 
find it to be negligible. Therefore our conclusion is based 
on the accuracy of CE04 clouds distribution, and is mainly 
limited by the XRT resolution and sensitivity. 

In the second (Section 4.2) we relax our assumptions re- 
garding the distance to the scattering clouds, and motivated 
by observations we assume a single predominant scattering 
screen. We keep the location of this screen a free parameter, 
and use Equation 5 to put simultaneous constraints on both 
the dust screen location and the dust properties, using a full 
data set including both spectral and temporal flux evolu- 
tion. The results of this analysis are limited mainly by the 
delay of the first observation and by the sensitivity of the 
XRT detector. 

4.1 Distance estimate based on known molecular 
clouds along the line of sight 

Equation 2 expresses the observed scattered flux as a func- 
tion of: (1) The scatterer's location relative to the source 
(i.e., x) and (2) Properties of the scattering dust. Therefore, 
if the location of the scatterers and the properties of the 
dust are known, one can extract the distance to the source. 
In fact, as we show here, it is enough to know the location of 
the dust screens along the line of sight and their dust column 
densities to put strong constraints on the source distance, 
even if the grain size distribution, Na is not well known. 
CE04 provide the locations of the scattering clouds along 



the line of sight to SGR 1806—20, as well as their optical 
extinction Ay, which indicate on their dust column densi- 
ties. Using this data we calculate the scattered flux radial 
proflle for different locations of the source. We then compare 
these "synthetic" profiles with the observed profiles and use 
the best fit profile to infer the distance. 

Table 2 lists properties of the molecular clouds along 
the line of sight, adapted from Table 1 of CE04. Following 
radial velocity measurements, each cloud has two distance 
solutions and was therefore attributed with a near distance 
and a far distance, as well as optical extinction Av- For 
eight clouds CE04 ruled out one of these distances due to 
additional considerations. We assume that the Ay given for 
each cloud is proportional to its share of the clouds' total 
dust column density (i.e. Ay oc A in Equation 6). Equipped 
with the clouds' positions and dust shares, we assume a 
single dust grain size distribution across all clouds and set 
cimax ~ 0.3 /im, amin = /im, q — —3. Assuming a source 
distance we use Equations 1 and 5 to calculate the scat- 
tered flux from each of the intervening clouds and construct 
a synthetic flux radial profile at any given time and energy 
band. In order to compare the profile calculated using Equa- 
tions 1 and 5 to the observations, it must be smeared by the 
PSF, the exposure duration and the width of the scattering 
clouds. We therefore introduce the Swift /XI{T PSF func- 
tion (Moretti et al. 2005), and we account for the observa- 
tions actual exposure durations to widen the profile beyond 
the PSF effect. Assuming a typical molecular cloud width 
of ~ 100 pc, the consequent smear is negligible. We repeat 
this calculation for an SGR 1806—20 distance that vary in 
the range 5—25 kpc with a 0.1 kpc resolution. Note that the 
synthetic profiles are calculated up to a normalization factor 
per energy band, accounting for the unknown direct fiuence 
from the source Se, as well as both Cdust and the constant 
of proportionality in ^ oc Ay. 

The actual observed signal contains two additional fea- 
tures beyond the dust scattered fiux, and we therefore model 
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the observed radial flux profile as a sum of the three fol- 
lowing contributions: (a) the synthetic profile of the dust 
scattered signal, calculated as described above (b) the mea- 
sured background profile, and (c) the halo around the SGR 
described in Section 2. Therefore, our model contains the fol- 
lowing free parameters: the halo amplitude and power-law 
per observed profile, the normalization factor of the scat- 
tered signal per energy band, and the SGR distance. 

We compare this model with the observed profiles from 
the first two epochs in Table 1, where the observed dust 
scattered ring is most evident, with each epoch split into 
two energy bands. For each distance we use a fit to find 
the halo amplitude and power-law per observed profile, and 
the normalization factor per energy band. We then calculate 
as a function of distance and choose the distance yielding 
the overall minimal the best fit distance. 

Three of the clouds lack an unambiguous distance de- 
termination and have instead both a near and a far distance 
estimates, similar for all three clouds. To account for this 
uncertainty we prepare two sets of scattered signal profiles 
as a function of SGR distance, one with all three clouds set 
at their 'near' position, the other with all three at their 'far' 
position, and repeat the above procedure for both sets. We 
verify that other combinations, where some of the clouds are 
at an opposite location, result in intermediate SGR locations 
and therefore the 'far' and 'near' configurations bracket the 
constraints on the location of SGR 1806—20. 

In order to estimate the distance error we use Monte 
Carlo simulations (following the prescription in Press et al. 
1986). For each Monte Carlo realization we generate a 
pseudo binned radial profile by drawing the "observed" ra- 
dial bins photon counts from a Poisson distribution around 
the best fit radial profile. We further prepare synthetic scat- 
tered profiles per SGR distance by drawing clouds' posi- 
tions^ and optical extinctions (from Table 2) assuming an 
approximately normal distribution for both. We repeat the 
fitting procedure described above for each of 10,000 Monte 
Carlo realizations and thus obtain the probability distribu- 
tion of the SGR distance (see Figure 2). 

To account for the dependence of the distance on the 
dust grain size distribution we run the complete procedure 
described above for nine different parameter choices of Equa- 
tion 5, with the permutations of Omax ~ 0.15, 0.3, 0.6 ^m and 
q = —3.5, —3, —2.5. We find the best fit distance scatter due 
to the model choice to be smaller than the scatter due to 
the errors of the binned signal and of the clouds' positions 
and opacities (a scatter of ~ 1 kpc between the two extreme 
model parameter choices, compared to ~ 3.5 kpc (la) scat- 
ter due to the errors). We therefore set the intermediate 
values of flmax = 0.3 fim and 5 = — 3 as our model parame- 
ters for the analysis. In Section 4.2 we measure q and find 
a lower limit for Omax, which are consistent with the pa- 
rameters we choose for this analysis, and which errors are 



^ CE04 use a systematic error of itlOkm s~^ on their measured 
velocity of the clouds in order to derive > 2cr distance errors. 
We therefore adopt a ±5km la velocity error for estimating 
the distance error to all the clouds using the Galactic rotation 
curve toward the line of sight presented in Figure 5 of Figer et al. 
(2004), based on Brand & Blitz (1993). The resulting distance 
errors are presented in Table 2. 
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Figure 2. The distributions of the distance to SGR 1806-20 
obtained by Monte Carlo simulations, with 68% and 90% con- 
fidence limits marked. The left (right) distribution matches the 
'near' ('far') clouds configuration, where the three clouds with am- 
biguous locations are fixed at their near (far) locations (see Table 
2). The best fit distance is IG^g'g kpc for the 'near' configuration 
and 11.9^2 5 '^P'^ the 'far' configuration (90% confidence). 



consistent with the range we use for this parameter inde- 
pendency check. The absence of a measured upper limit for 
cimax has no infiuence on this analysis since the dominant 
clouds fall within the power-law regime (see Equation 4) for 
any flmax bigger than the value chosen here. 

The effect of the SGR position on the scattered sig- 
nal and on the fit quality is demonstrated in Figure 3. 
For an SGR at 16 kpc, the clouds concentrations at 4.5 
and 6 kpc nearly merge, overlapping with the measured 
signal, while the 3 kpc clouds are negligible. In contrast, 
with the SGR at 6 kpc the clouds concentrations are more 
separated, with the dominant 4.5 kpc cloud and the mea- 
sured signal largely unsynchronized. For the 'near' con- 
figuration, the best fit distance is 16l3 g '^P^ (90% con- 
fidence), with x^/dof = 43.3/47. For the 'far' configura- 
tion, best fit distance is 11.91^2 5 ^P^- (90% confidence) with 
X^/dof = 45.9/47. Since the goodness of fit for both con- 
figurations is similar we use their combination to set con- 
servative limits, with a lower limit of 9.4 kpc and an upper 
limits of 18.6 kpc on the distance to SGR 1806-20 at a 90% 
confidence level. 

We note that taking the main bulk of the dust along 
the line of sight according to CE04, concentrated at 4.5 kpc, 
and substituting in Equation 1 with our observed ring's d 
and t, results in a distance of 15 kpc for SGR 1806—20. Our 
above analysis wraps this crude estimation with the con- 
fidence limits permitted by the XRT observations and the 
CE04 errors. 

The weakest link of our analysis is the accuracy of 
the clouds' locations as reported in CE04. As an example, 
Figer et al. (2004) suggested, based on a different method 
for measuring the radial velocity, a distance of 11.8 kpc (vs. 
15.1±0.9 of CE04) to cloud MC13A. Although this cloud 
has minor influence on our analysis, the discrepancy in this 
measurement demonstrates the limitation of our method. 
However, the goodness of fit we obtain, x^/dof ~ 1, pro- 
vides a consistency check for the CE04 clouds distribution 
along the line of sight. In Section 5 we argue that future 
post-burst observations could resolve the location of a few 
individual clouds. This would potentially permit an SGR 
distance estimate that depends, independently, on the lo- 
cation of each of the resolved clouds, thereby resulting in 
a more robust and accurate distance measurement for this 
SGR. 



SGR 1806—20 distance and dust properties in molecular clouds 7 




Figure 3. A comparison between the best fit obtained for SGR 1806—20 located at 6 icpc vs. 16 kpc, witii clouds at the 'near' clouds 
configuration (see Table 2). Circles mark the best fit model binnned photon counts while dots with error bars (Poisson, Icr) mark the 
measured bin photon counts. Arrows indicate the positions in kpc of the main clouds concentrations corresponding to each peak. An 
SGR distance of 16 kpc evidently provides a good fit to the data, while a distance of 6 kpc is ruled out at a > 3(t level. Model components: 
Dotted line is background, dotted-dashed is the fading halo, dashed is the calculated signal profile, continuous line is sum of the three 
components of the best fit model. 



4.2 Constraints on dust properties 

Here we assume a single predominant dust screen at an un- 
known a-priori location, thus dropping our previous rely on 
CE04 results. This assumption is motivated by the fact that 
the XRT observations show a single ring. Fitting the ob- 
served radial profile fluxes and photon energies using Equa- 
tion 5, we obtain best fit values for flmax, q and B{E). Since 
fimax ~ flmax/ (1 — x) the fit rcsults in a joint constraint on 
flmax and X. The model we use assumes a power-law grain 
size distribution A'^a oc a' with grain size flmin < a < Omax, 
leading to the three scattered flux regimes described in 
Equation 8. Thus, constraining q requires observations at 
epochs or angles within the power-law regime, while con- 
straining flmax requires the further inclusion of earlier obser- 
vations, towards the constant flux asymptote. 

Fitting the observations to the dust model is done in 
two steps. First we find the integrated scattered signal flux 
within the ring and its associated flux probability distribu- 
tion (from simulations) per epoch (for the four epoches in 
Table 1) per energy band by fltting the observed data to a 
model. Then we use these fluxes as a function of epoch and 
energy, and their probability distributions, to obtain best 
flt values for fimax, q and B{E) by means of a Maximum 
likelihood flt using Equation 5. 

We model the dust scattered signal as a Gaussian pro- 
flle, thus accounting for Swift PSF and exposure duration 
effect, assuming that the effect of the clouds' width is negli- 
gible. In order to calculate the scattered signal net flux per 
energy band per epoch, we fit each observed radial profile 
to a model composed of the following three components: (a) 
the Gaussian profile representing the scattered signal (b) 
the measured background profile, and (c) the halo around 



the SGR described in Section 2. The contribution of each 
component to each bin is weighted according to the bin's 
average exposure time due to the exposure maps. We run a 
minimum fit with the fading halo amplitude and power, 
and the Gaussian amplitude and width as the free param- 
eters. For the first observation we also fit for the radius of 
the Gaussian center and obtain** 6 = 119 ± 5" {2a). We fix 
the ring location of the other observations using 6 [t) oc \/i 
(Equation 1). Figure 4 shows the best fit model for the first 
two epochs. The best fit Gaussian area translates to the pho- 
ton count per second related to the scattered signal, which 
we convert to fiux using NASA's HEASARC web based tool 
WebPIMMS^ 

We find probability distributions of each ring's inte- 
grated fiux by Monte Carlo simulations. The best fit ra- 
dial profile from the fit is used as the base data set for 
preparing (assuming Poisson distribution for the bins pho- 
ton count) 5000 binned profile realizations. Collecting the 
best fit fiux of each realization we compile the flux proba- 
bility distribution. 

These distributions are then used as the input for a 
Maximum Likelihood fit, where we flnd the best flt values 
for fimax and q, along with two normalization factors Be, 
one per energy band. For a given set of parameters ayn!ix,q 

* An upper limit of 6.37 kpc for the dominant dust screen, re- 
gardless of the source position, is obtained by substituting the 
observed ring's 9 and t in Equation 1 (a 6.6 kpc upper limit was 
reported by Goad et al. 2006). 

^ Using a power-law source model with photon index -1, i.e. a 
blackbody Rayleigh- Jeans law (e.g. Feroci et al. 2004), and Nh = 
6 X 10^2 cm-2 (Sonobe et al. 1994). 
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Figure 4. fits for the first two epochs in the two energy bands. 
The dotted line is background, dotted-dashed is the fading halo, 
bottom continuous line is the scattered signal Gaussian and top 
line is the sum of the three components. Dots with error bars 
(Poisson, 1(t) are the binned measured photon counts. The best 
fit Gaussian area translates to the flux per epoch and energy band 
(Table 3), flux errors distribution are found using Monte Carlo 
simulations (see text). 
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Figure 5. Maximum likelihood best fit model (Equation 5). 
Lower line and accompanying dots and error bars (Table 3) are 
the 2.5—4.5 keV band. The best fit result for the power-law index 
q is — 3.3^Q^ (l""); with errors determined from Monte Carlo sim- 
ulations. The best fit result for Smax = iniax/(l — is 0.33 fim. 
Since the first observation is almost completely within the flux 
power-law regime (see Equation 8), we can not set an upper limit 
for Smax- The 2cr lower limit is 0.21 fim (see Figure 6). 



and B{E) we use Equation 5 to calculate the flux at each 
observation epoch and energy band. We then evaluate the 
likelihood of this set using the measured flux distribution. 
The set of parameter values which yields the highest product 
of likelihoods for all epochs and energy bands is the best fit. 
Figure 5 shows the Icr bars of the flux distributions, and the 
corresponding best fit model. 

Estimation of the confidence limits of these parameters 
is calculated by Monte Carlo simulations. For each Monte 
Carlo run we generate a realization of the fiuxes in Table 3, 
with the ring's integrated flux values drawn from the actual 
best fit distribution in Table 3. The errors around the drawn 
fiux values are taken to be the same as the original errors in 
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q 

Figure 6. 68% and 90% likelihood contours of q and Omax as 
determined from the Monte Carlo simulations. Note the lack of 
an upper limit for Smax- The marginalized distribution of each 
parameter are shown above (q) and to the right (Smax). The lines 
mark the 68% confidence limits for each single parameter, based 
on an integration of the probabilities for Smax ^ 3/im. 

Table 3. We fit 5000 Monte Carlo realizations for fimax, q and 
B{E), thus obtaining the probability distribution of these 
parameters. 

Fluxes used for the maximum likelihood fit, with la 
errors and values, are listed in Table 1. The best fit result 
for q is —3.3^0 7 consistent with the commonly cited 

—3.5 from Mathis et al. (1977) and with an implicit —3 of 
Draine (2003). The best fit result for Smax is 0.33 /^m. As 
indicated by the likelihood contours of q and Omax shown in 
Figure 6, we are unable to constrain Omax from above, which 
is the result of the earliest observation being too adjacent 
to the power-law asymptote regime (see Figure 5). We thus 
find only a lower limit of 0.21 fim {2a) for Omax- 

In our fit we constrain Smax, which is a combination of 
cimax and X (ciinax = ° ) . lu Order to learn further about 
either x or Omax one needs to have prior knowledge of ei- 
ther ttmax or X, respectively. Since 6 (t) is known, a mea- 
surement of X will provide an SGR distance estimate that 
is independent of our previous analysis. However, prior con- 
straints on the distance to SGR 1806—20, 6—15 kpc (Sec- 
tion 1), translate, together with the locations of the ring, 
to 0.30 < X < 0.51 (Equation 1). Therefore, improving 
these constraints on x requires constraints on both Smax 
and Omax that are better than 25%. The uncertainty in the 
value of Smax in the diffuse ISM is larger, while its value 
in molecular clouds is practically unknown. Therefore, our 
measured amax, and even an ideal errorless fimax, are use- 
less for improving the constraints on the distance to SGR 
1806—20. This conclusion is similar to that of Tiengo et al. 
(2010), which find that the distance to the X-ray pulsar 
1E1547. 0-5408 cannot be conclusively determined due to the 
uncertainty in the dust properties, although they obtained 
superb data. 

In the case of SGR 1806-20, however, there are robust 
(although rather loose) limits on its distance. Thus, the prior 
constraints on x can be used to efficiently constrain Omax, 
which is of special interest in our case of molecular clouds. 
The robust constraint of d > 6 kpc for the SGR's distance 
(Cameron et al. 2005) matches to a; < 0.51, which translates 
to flmax > 0.1 ^m {2a). 
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Table 3. Best fit ring's integrated fluxes used for the maximum likelihood flt 
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Epoch 


Time since burst 


Energy band 


Flux 




xVdof 




[hr] 


[keV] 


[ergcm~^s^^] 






1 


30.8 


2.5-4.5 


1.63 ±0.27 X 10" 


-12 


2.95/4 


1 


30.8 


4.5-8.5 


4.45 ± 1.13 X 10" 


-12 


2.00/5 


2 


93.5 


2.5-4.5 


3.39 ±0.71 X 10" 


-13 


6.15/15 


2 


93.5 


4.5-8.5 


5.73 ± 2.45 X 10" 


-13 


17.25/15 


3 


121.1 


2.5-4.5 


2.88 ±8.70 X 10- 


-14 


7.59/11 


3 


121.1 


4.5-8.5 


4.43 ±2.92 X 10" 


-13 


13.31/11 


4 


147.6 


4.5-8.5 


3.64 ±3.12 X 10" 


-13 


16.83/11 



The best fit results for the normahzation factors are 
B(2.5 - 4.5keV) = 2.89lEJ;|^ ^ ^^'^^ ergcm"^ s"^ asec"^ 
B(4.5 - 8.5keV) = x 10"^^ ergcm"^ s"^ asec"^ 

As a consistency check, we use the ratio between the two 
normalization constants to find the Hydrogen column den- 
sity towards SGR 1806—20. The normalization constant de- 
pends on photon energy through Se, the direct observed 
fluence. Given the source spectrum and the x-ray extinction 
cross section per H nucleon per energy ct{E), Nh can be 
extracted using 5*^ oc exp{—a{E)NH}- Assuming the x-ray 
extinction is mainly dust driven, we take the x-ray extinc- 
tion cross section per H nucleon due to dust for each energy 
band from Draine (2003) (assuming Rv ~ 3.1). For the burst 
x-ray spectrum at our energy bands we follow Feroci et al. 
(2004) and Olive et al. (2003) in adopting a double black 
body spectrum instead of the OTTB spectrum commonly 
used to describe SGR bursts in the range 20 — 200 keV 
(e.g. Golenetskii et al. 2005). Using the average double black 
body spectrum of kTi = 3.4, kT2 = 9.3, L1/L2 = 0.85 
from Feroci et al. (2004) with our best fit ratio we get 
Nh = 9.64j;4^ X lO^^cm"^, consistent with Sonobe et al. 
(1994). 

Reversing the last procedure we find the normalizations 
ratio between the 2.5—4.5 keV band and the 0.5—2.5 keV 
band, which we did not use due to lack of signal, to be 
1000:1, explaining the lack of signal beyond the background 
noise in the 0.5—2.5 keV energy band. 



5 FUTURE SGRS OBSERVATIONS 

The accuracy of the derived distance to SGR 1806—20 and 
the properties of the dust along the line of sight are re- 
stricted by large error bars due to the weak signal and low 
resolution of the Swift fKKT , and by the long delay between 
the burst and the first echo observation. Nevertheless, our 
analyses demonstrate the wealth of information that can 
be extracted from observing dust scattered x-ray echoes of 
bursts. Encouraged by these results, we provide guidelines 
for optimal future observations of magnetar bursts, showing 
that improved observations' timing, resolution and sensitiv- 
ity would yield far better constraints on both the distance 
to an SGR and the properties of the intervening dust. 

In order to evaluate the expected angular scattered 
width of the ring we compare the effects of clouds width, 
exposure duration and PSF. Substituting an SGR distance 
of 15 kpc and x = 0.3 (i.e. main cloud at 4.5 kpc) in Equation 
1 we get a spread of roughly 1" for an assumed typical cloud 



width of 100 pc, and a spread of AO = 0.36" \^t^\^,t ^ iLdj 

due to exposure duration (where t is the time passed in sec- 
onds since the burst direct observation). This yields a width 
of ~l" for the first XRT epoch, at a delay of 30 hours and 
duration of 40 minutes, due to the non-PSF effects. Thus 
the width is practically set by the PSF of the telescope, 
with ~20" for Smft /XRT vs. 1" for Chandra. 

Effective constraining of amax requires observations ear- 
lier than 30 hours after the burst, and such a response time 
is suitable for Sm/t/XRT. With a stronger scattered sig- 
nal and a sharper decay of the King component for earlier 
observations, resolving the scattered signal from the source 
should be possible for a ~1 hr observation taken 12 hours 
after the burst (^obs = 75"). Since observations should also 
cover the scattered flux power-law regime, the burst fluence 
should be similar or larger than the 2006 August 6 burst. 

Chandra observation can be obtained with a response 
time of a few days. The improved resolution and collection 
area of Chandra could significantly reduce errors and im- 
prove statistics. Chandra also enables a continuous observa- 
tion due to its high orbit, thus minimizing the spread due 
to exposure duration. A resolution of 2"— 3" (accounting 
for PSF, exposure duration and cloud width) and a signal 
roughly four times that of Swift should improve the signal 
to noise ratio, and enable resolving the location of individ- 
ual clouds along the line of sight. Each of the resolved clouds 
will provide an independent constraint on the SGR distance. 
This will provide a more accurate and, more importantly, 
a more robust distance estimate. In Figure 7 we illustrate 
the advantages of using Chandra, comparing what Chandra 
would have seen at the second XRT epoch to the actual 
Swift observation. 

The burst energy determines the proper time for the 
last effective observation. A delay of 30 hours, equivalent 
to 9 = 120", is located within the power-law regime fiux 
decay (Equation 8), as seen in Figure 5. Using our q — —3.3, 
combined with our measured scattered fiux, we get 



-Fscat — Fq 



/ ^ delay \ 
Augusts, 2006 V^O^^/ 



(10) 



1.63 X 10"^^ergcm"2s"^ and 
ergcm"^ s~^ are the fiuxes 



where Fo(2.5 - 4.5 keV) 
Fo(2.5 - 4.5 keV) = 4.45 x 10"^^ 
from table 3. The equation is valid for tdoiay J5 30 hours 
(power-law regime). We calibrate the background photon 
count by the total background photon count measured for 
Swift (e.g. 0.14 phs~^ for 2.5—4.5 keV), with the background 
per radial bin piling up as oc By requiring the bin 
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which contains the largest fraction / ^ 1 of the scattered 
signal to be times larger than the background error we 
get a rough estimation of the maximum delay for an effective 
observation: 

t = 

, / 5 . / 10" A. T^r- ATA''-^^ (11) 

\ ^ August 6, 2006 NtT-' y bin width A^wift i hr ^ 

where to(2.5-4.5 keV) = 135 brand fo(4.5-8.5keV) = 
120 hr, Acf! is the detector's effective area for the chosen 
energy band, and At is the observation's exposure time in 
hours. 

For the Swift case at 2.5 — 4.5 keV we get a > 5(t de- 
tection for the first observation, at t=30hr, and a ~ 5a 
detection for the t=100hr second observation. A Chandra 
■A-cB ~ 4j4^g'^', At = 5hr (for the latest observation), a bin 
width of 3" and / ~ 1 permits a 5cr detection as long as 
370 hr after a burst of the same fiuence, thus probing dust 
grains down to a size of ~ 0.022 fim (vs. ~ 0.035 /im in our 
case). Alternatively, a Chandra observation of At = 3hr af- 
ter 100 hr can yield a 5cr detection for a burst of a fiuence 
ten times weaker than that of 2006 August 6. Such SGR 
bursts are more common, offering more opportunities for 
observations. 

Equations 10 and 11 quote the fiuence of the 2006 
August 6 burst. The measured event fiuence at the range 
20-200 keV is $ = 2.4 x lO'^'erg • cm'^ (Golenetskii et al. 
2006). We assume, at least for minor/intermediate bursts, 
a similar spectrum across different energy outputs (e.g. 
Olive et al. 2003; Feroci et al. 2004). Hence, for events which 
fiuence is given at 20—200 keV, one should substitute the 
above cited fiuence along with the new event fiuence. For 
events where the band fiuence (e.g. 2.5 — 4.5 keV) is more 
accessible, or events for which the similar spectrum assump- 
tion does not necessarily hold (e.g. giant flares), we sup- 
ply an estimate of the 2006 August 6 burst fiuence at the 
2.5 — 4.5 keV and the 4.5 — 8.5 keV energy bands. Olive et al. 
(2004), Feroci et al. (2004), Nakagawa et al. (2007) and 
Israel et al. (2008) obtain satisfying fits for the spectra of 
minor/intermediate SGR bursts using models composed of 
two blackbody components, a harder kT « 9 — 10 one 
and a softer kT « 3 — 4 one. Israel et al. (2008) compared 
three of the cited works and found the bolometric lumi- 
nosities of the two components to be similar, with a pos- 
sible saturation of the soft component above lO'^^ergs"^. 
We adopt an equal luminosity for both components. Cali- 
brating all the cited models with the SGR fiuence per keV 
at kT = 20 (derived from the total measured fiuence and an 
OTTB spectrum of kT = 20 keV over the range of 20-200 
keV), we calculate the band fiuence according to the black- 
body temperatures of each cited model, and estimate the 
fiuence as the geometric mean of the two extreme model 
results, obtaining $2.5-4.5 = LO^q I ^ lO"'^ erg cm~'^ and 
$4.5-8.5 = 8.313 2 X 10~^ ergcm"'^. These values can be used 
for equations 10 and 11 along with the respective band fiu- 
ence of a new event, regardless of the event's spectrum. 

A burst brighter than the one we analyzed should en- 
able probing smaller dust grains and even, if bright enough, 
allow the determination of amin, the lower cutoff of the dust 
grain size distribution. The launch of NuSTAR next year 
may provide an opportunity to constrain a^iin by observing 
echoes from an intermediate burst, since for a given amin 
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Figure 7. A comparison of the best fit photon counts for the 
second epoch 2.5— 4.5keV band with SGR at 16kpc, Swift vs. 
Chandra. Based on a Chandra 4-fold signal and 3" resolution. 
Contributions of clouds scattered signal, background and fading 
halo are shown as continuous lines. Poisson la confidence limits 
for the best fit photon counts are marked as gray filled area. Swift 
actual measured photon counts are marked in circles. 

the transition time from the power-law regime to the ex- 
ponential decay regime is proportional to 6^ and therefore 
to (Equation 4). NuSTAR has a sensitivity range of 

5 — 80 keV and a larger effective area compared to Chandra 
(for 15 < 35 keV ) . In addition, the signal to noise at higher 
energies should improve due to the spectrum of SGR bursts 
and to the decrease in the background level. On the other 
hand, NuSTAR's relatively low resolution of 43" will dilute 
the signal. Taking it all into account, NuSTAR should be 
able to follow hard x-ray echoes from intermediate bursts 
for a few days after the burst, thus probing much siualler 
grains, and potentially identifying flmin. 

Generally, x-ray echoes observations taken during dif- 
ferent bursts of the same SGR can be combined to improve 
signal to noise ratio and thus improve the constraints dis- 
cussed in this work. 



6 CONCLUSIONS 

We used Swift /XI{T observations following an intermedi- 
ate burst of SGR 1806—20 to constrain the distance to the 
SGR and to find the first x-ray echo constraints on molecu- 
lar clouds' dust properties. Based on the molecular clouds' 
properties along the line of sight found by GE04 and the 
observed echoes of the burst we constrain the distance with 
a lower limit of 9.4 kpc and an upper limits of 18.6 kpc at 
a 90% confidence level. This upper limit can be considered 
the first direct upper limit set for the distance to the SGR, 
as the upper limit set by Cameron et al. (2005) is question- 
able, and other distance estimates are based on association 
rather than measured emission from the source. Our dis- 
tance constraints favor an energy output of ~ 10*^ergs~^ 
for the 2004 giant flare of SGR 1806-20, leaving Galactic vs. 
extragalactic giant bursts rates possible tension and the lack 
of post-burst emission features changes in SGR 1806—20 as 
open issues. 

We introduce the use of observations of dust x-ray 
echoes for probing the dust properties of molecular clouds. 
Fitting the spectral and temporal signal evolution using a 
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dust scattering model with an assumed power-law dust grain 
size distribution we find a power-law index of g = — 3.31q'7 
(la) and that the dust maximal grain size amax > 0.1pm 
(2cr). These results are of special interest since the dust along 
the line of sight to SGR 1806—20 resides in molecular clouds, 
which dust properties have been poorly explored. The con- 
straints we got, imply that the dust grain size distribution in 
molecular clouds may be similar to the one found in diffused 
ISM (e.g. Mathis et al. 1977). 

The wealth of data obtained using Suii/i/XRT encour- 
age us to suggest the use of future x-ray echoes observations. 
For an intermediate SGR 1806—20 burst, a Swift observation 
starting half a day after the burst and a Chandra observa- 
tion a few days after the burst should yield superior dust 
grain size distribution characterization and SGR distance 
constraints, respectively. A Chandra observation should en- 
able, in addition, a quality mapping of clouds' locations 
along the line of sight. NuSTAR, planned to be launched 
next year, could potentially probe Omin, the minimal dust 
grain size. 
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